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type system are investigated. This study is performed from the point of view 
of a two color dressed state basis. Both of these processes are possible due to 
atomic coherence and quantum interference contrary to simple intuitive predic- 



' tions. Merely on physical basis, one would expect a complementary process to 

the amplification without inversion. We believe that absorption in the presence 

^ of population inversion found in the dressed state picture utilized in this study, 

00 ' 

^ | constitutes such a process. 

Q_i| Novel approximate analytic time dependent solutions, for coherences and 

+-> ' 

populations are obtained, and are compared with full numerical solutions ex- 

d ' 

qh hibiting excellent agreement. Steady state quantities are also calculated, and 



the conditions under which the system exhibits absorption and gain with and 



without inversion, in the dressed state representation are derived. It is found 
that for a weak probe laser field absorption with inversion and amplification 
without inversion may occur, for a range of system parameters. 
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I. INTRODUCTION 



Recently there has been tremendous interest in the study of light amplification and lasing 
without the requirement of population inversion (LWI), potentially capable of extending the 
range of laser devices to spectral region in which, for various reasons, population inversion is 
difficult to achieve. These spectral regions include UV which can be obtained from atomic 
vapor, and mid-to-far infrared, obtained by interssubband transitions in quantum wells. Many 
models for LWI have been proposed, mostly three and four level schemes, in A and V configura- 
tions. The dependence of optical gain on system parameters has also been investigated PHH]]. 
The key mechanism, which is common to most of the proposed schemes, is the utilization of 
external coherent fields, that induce quantum coherence and interference in multilevel systems. 
An exception of LWI without the use of coherent fields was also reported ||. In particular it 
was shown that if atomic coherence between certain atomic states is established , different ab- 
sorption processes may interfere destructively leading to the reduction or even the cancellation 
of absorption ]7,13]. At the same time stimulated emission may remain intact, leading to the 
possibility of gain even if only a small fraction of the population is in the excited state. 
Experimental observations of inversionless gain and lasing without inversion have been reported 
by several groups |T8|-F2T|. Inversionless lasers have been shown to have unique properties such 
as non classical photon statistics and substantially narrow spectral features fl2"2] , |2"5|] . In a recent 
paper, Y. Zhu |12[ analyzed the transient and steady state properties of light amplification 
without population inversion in a closed three level V type system in the bare state basis. 
Steady state dressed state populations were also calculated, in the limit of strong driving laser. 



In this paper we study a V- type three level model within the framework of the dressed state 
basis, and give explicit analytic time dependent solutions, as well as steady state solutions for 
populations and coherences. This paper details the calculations of LWI in the dressed state 
basis, which is valid to some extent for atoms dressed both by the pump and probe lasers. 
More over we offer a set of theoretical tools allowing one to obtain explicit time dependent 
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solutions for populations and coherences. These calculations show the possibility of inversion 
and inversionless gain in the dressed basis as well as gain without population inversion in 
the bare state representation. A novel interesting feature found in this study is the existence 
of absorption despite population inversion. Although this effect is contrary to simple physical 
intuitive explanation of absorption, this process has a conceptual reasoning. From basic physical 
arguments one expects a complementary process to amplification without inversion. We believe 
that absorption in the presence of population inversion, found in the dressed state picture, 
constitutes such a process. It is another manifestation of the quantum interference that may 
occur in multilevel systems where coherently prepared states present the possibility of interfering 



channels [p4 |. This phenomenon should be interpreted as a quantum interference constructive 
process for the stimulated absorption just as LWI is obtained from a quantum interference 
destructive process of absorption. 

In Sec. II we present the model system and the master equation used to derive the equations 
of motion for the elements of the density matrix. We have chosen to employ a fully quantum 
mechanical Hamiltonian, even though later on the density matrix equations are reduced to 
their semiclassical version. The quantum mechanical Hamiltonian gives rise to a simple picture 
of the stationary dressed states. In Sec. Ill the dressed levels are introduced and the master 
equation is projected over the dressed state basis. Physical interpretation of relaxation coeffi- 
cients in the dressed basis is given and the role that quantum coherences and interferences play 
is elucidated. In Sec. IV we present time dependent and steady state solutions for the dressed 
state populations and coherences. Comparison is made with the bare state results. 
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II. HAMILTONIAN AND MASTER EQUATION FOR THE ATOM 



Let us consider the closed V-type three level system illustrated in figure III The transition 
| a >^-> \b > of frequency Ub a is driven by a strong single mode laser of frequency u>l. The tran- 
sition \a >^ \c > of frequency uo ca is pumped incoherently with a rate A. A single mode probe 
laser (not necessarily weak) is applied to the transition \a >^-> |c >. 7& ( r y c ) is the spontaneous 
emission rate from the state \b > (|c >). The states \b > and |c > are not directly coupled. 



We have chosen to work within the frame of the master equation for the atom, since it being 
an operator equation independent of representation, it can be projected over any basis. We use 
a generalization of a standard master equation adjusted to account for the scheme described 



above f25fl . The master equation is given by: 

o = - ^[-^> °] ~~ ~2 l s + s - a + o"s + s_] + 7 6 s_crs + - —[s' + s'_(T + as' + s'_) + ^ c s'_as' + 

— — [s' + s'_a + as' + s'_] — — [s'_s' + a + us'_s' + ] + As' + as'_ + As'_as' + , (1) 



here a is the density operator for the atom. s + (s_), s' + (s'_) are the atomic raising and 
lowering operators, for the \a >^ \b > and \a \c > transitions respectively. H is the 
Hamiltonian of the global system and we take it to be fully quantum mechanical. The quantum 
mechanical Hamiltonian gives rise to a simple picture of the stationary dressed states. The 
Hamiltonian in the dipole and rotating wave approximation is given by: 

H = hu ba \b >< b\ + hui ca \c >< c\ + huj L (a)a + -) + hu p (a^a' + -) 

+ g(s + a + s_a j ) + g'(s' + a' + s'_a' ] ). (2) 

g and g' are coupling constants and are assumed to be real. The eigenstates of the unperturbed 
part of the Hamiltonian form a three dimensional manifold labeled by the atomic indexes, the 
laser photon number N and by the probe photon number N'. The manifold is written 

e(N, N') = {\a, N + 1, N' + 1 >, |6, N, N' + 1 >, |c, N + 1,N' >}. (3) 
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We represent the uncoupled eigenstates of the atom and the two non-interacting field modes 



as: 



\a,N + l,N' + 1 >-- 



v / 



\b,N,N' + 1 >= 



V / 



\c,N + l,N' >= 



V 1 / 



(4) 



In this basis the hamiltonian takes the form 

/ 



H 



-Q -G 

In, -Q -Al • (5) 

G -A 2 

Where we have defined the Rabi frequencies and the detunigs, in their quantum form by 



V 



/ 



gVN + 1 = hVL ; —g'\/N' + 1 = KG ; A 1 = uj l - uo ba ; A 2 = uo p - uo ca . 



(6) 



To obtain the semiclassical equations of motion for the elements of the density matrix we 
project the master equation (HD over the basis @ and perform the following reduction operation 
by introducing the reduced populations and coherences pij via: 

Paa= E < a, N + 1, N' + l\a\a, N + 1,N' + 1>, (7a) 

(7b) 



N,N> 



p ab = E < a, N + 1, N' + l\a\b, N, N' + 1 > . 

N,N> 

and similar relations for the other populations and coherences. Taking into account the above 
reduced quantities we obtain the density matrix equations 



Paa = -Apaa + (A + 7 c )p cc + %p bb + itt(p ba ~ Pab) + iG(p ca ~ p ac ), 
Pbb = -IbPbb + i^iPab ~ Pba), 



Pc 



Ap aa - (A + 7 c )p cc + iG(p ac - p ca ), 



Pab = -[-(A + 7 fe ) + iA^Pab + iQ(p bb - paa) + iGp cb , 
Pac = _ [^(A + 7c) + «A 2 ]p ac + lG(p cc - paa) + i^Pbc, 

pbc = ~[^(A + lb + 7c) + ^(A 2 - Ai)]p 6c + iVtpac - iGp ba . 



(8b) 



(8c) 



(8d) 



(8f) 



III. DRESSED STATES AND DENSITY MATRIX EQUATIONS IN THE DRESSED 

STATE BASIS 



The dressed states are obtained by finding the the eigenvectors of the interaction Hamiltonian, 
Eq. (|5p. To simplify things a little we take both the driving laser field and the probe field to be in 
exact resonance with the corresponding bare state transitions, i.e., we take A x = A 2 = 0. When 
the detunings are set to zero we notice that the energies in the bare state basis are all degenerate, 
and in fact equal to zero (in the interaction picture). Carrying out the diagonalization procedure 
we obtain the following eigenstates: 

\a(N)(N') > = -- \b, N,N' + 1>+^ \c, N+ 1, N' >, (9a) 
\P(N)(N')> =-^\a,N + l,N' + l>+^L\b,N,N' + l>+-^L\c,N+l,N'>, (9b) 

\-y(N)(N') > =^=\a,N + l,N' + l> 1 6, N, N' + 1 > +-J- \c, N + 1, N' > . (9c) 

and the corresponding energies 

-E| 7 > = — TiR ; E\ a> = ; -^|/3> = hR- (10) 

Where we have introduced the on resonance, two field generalized Rabi flopping frequency 



R = W + G 2 . 

Note that for the case G = equations ([T5|) and (^c|) render the usual coupling and non 
coupling dressed states, while the \a > state which is identical to the |c > state is not involved 
in the interaction altogether. The energy ladder is shown schematically in figure 2. We can see 
that one state remained intact, while the other two states were displaced by an energy amount 
of HR, with respect to the bare states. In the strong driving field limit, i.e., when Q » G, we 
see that the |a(A^)(A^') > state has the character of the excited bare state \c, N + 1, N' > and 
hence is expected to be less populated than the other two states. By contrast, the \P(N)(N') > 
and |7(iV)(iV') > states have a ground state character, and hence will be more populated than 
the \a(N)(N') > state. However both \/3 > and |7 > states are also contaminated by the same 



amount of the first excited level and thus they are expected to possess the same population 
content. 

The eigenstates of equations (pa|)-(p^) define a rotation matrix (transformation matrix) 



T 







g n 

' R R 



\ 



i n g 



V2 V2R V2R 



i n g 

V \/2 \/2R V2R J 



that diagonalizes the Hamiltonian of Eq.(§) via the matrix product THT . Thus the 
density operator in the dressed atom basis, p Dr , will be given by the matrix product 



P 



Dr 



Tp B T -l 



(12) 



where p B is the density operator in the bare basis. 

Projection of the master equation over the dressed state basis yields particularly simple 
equations for the first part of Eq. (|TJ) , i.e., the Hamiltonian part of the master equation. However, 
in the dressed atom basis the relaxation terms of Eq.(|Tj) give rise to equations that are not as 
simple as equations (|8a] -|3f|). In particular couplings between dressed state populations and 
coherences between two dressed states appear. In the next section we present an approximate 
version of the complete set of equations given bellow. 

The equation of motion for the density matrix elements in the dressed state representation 
are given by: 

Paa = -(T a + A')p aa + T(p a /3 + Pj3a) + f (Pa 7 + P 7 a) + ~A' ' ppp 
PaP = -(Tap - iR)PaP - (I> - -A')Pq7 - T'p/37 + (f - f')p~,p 



13a) 



+ TPaa + (f + T')ppp + f'p 77 , 



(13b) 



p ay = -(T a p + iR)por, - (Tp - -A')p a p + (f - T')pp y - f'p^p 



+ Tp aa + f'ppp + (f + P)p 77 (13c) 

ha = -( F /3 + ^ A )P/3/3 + ^(P* + A')Paa + (I> + - ~A')p 77 

- f(p Q7 + p 7Q ) + ~A'(pp 7 + p 7 p), (13d) 
P/3 7 = -(I> 7 + 2iR)pfr - (I> + ^ A - ^A')p 7/3 + f (p aP + p 7a ) 

+ 2f (pp a + p a7 ) - ~(r a + A')p aa - (2I> - -A')(ppp + p 77 ) (13e) 
p 77 = -(Tp + ^A)p 77 - T{p aP + p^) + -A'(p^ 7 + p 7/3 ) 

+ I(T a + A') PaQ + (I> + - IaOp/j/j (13f) 

Where we have again made use of the reduction operation 

PiJ = E < i{N){N')\*\j(N)(N') > i,j = \a >, \(3 >, | 7 > (14) 

N.N' 



In obtaining equations (|13a ) - ( |131j) we have introduced the following notation 



r> = r> + i(r„ + a + a'/2) , r> 7 = 3i> + ~a - a', 

* Gft , n - GQA . 1 AQ 2 

r = ^l^hb - 7 c - A) , P = — =— , A' 



r a , are the spontaneous emission decay rate of the \a >, \f3 > and | 7 > states (| 7 > 
also decays with a rate Tp). More precisely the state \a(N)(N') > decays by spontaneous 
emission with a rate T a to the levels \{3(N - 1)(N') >, \j{N - 1)(N') > , \0(N)(N' - 1) > 
and | 7 (iV)(A^ / - 1) >. Similarly the levels \/3(N)(N') >, \j(N)(N r ) > decay with the same rate 
F/3 to the same levels as \a(N)(N') >. The coherences p a p, and p ai (pp 7 ) decay with a rate 
T a p (r^). A' is a dressed picture pump rate which causes population and depopulation of the 
dressed levels. It also has an important influence on the coherences as can readily be seen from 
Eqs. (|13a|) - ( |13f|) . T and P are identified as interference terms. They involve the product of two 
Rabi frequencies. Both parameters vanish whenever either G or Q are zero. These terms are 
responsible for the amplification without inversion and for the absorption despite the inversion. 
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This fact is verified numerically. When we have set both T and f ' zero (this happens when 
A = and 7^ = j c ) any previously obtained gain has vanished. The first and second terms in 
T aj 3 describe damping of the atomic coherence due to radiative transitions of the levels involved 
to lower levels, and is equal to half the sum of all transition rates starting from \a(N)(N') > 
and \/3(N)(N') >. The third term in T a p describes coherence damping due to the incoherent 
pump. The interpretation of is similar except that the 31^ is composed of a 2Tp term 
responsible for the coherence damping due to radiative transition, plus a single Tp component 
resulting from transfer of coherence from higher levels belonging to higher manifolds. This fact 
would have been transparent had we written the non reduced version of Eqs. ( |13a|) - ( |13i|) . 



Inspection of equations (13) reveals that p a p and p ai have the same free frequency R, however 
they oscillate out of phase. The free evolution frequency of p^ is twice as large, as both levels 
\(3 > and I7 > are contaminated by the bare ground state \a >. Note that the closure of the 
system is satisfied by Eq. (|13a| )-(|131| ) , i.e., ^(p aa + Ppp + P 77 ) = 0. Gain or absorption coefficient 
for the \j >—>■ \i > transition is proportional to Im\ftij]. In our notation amplification will occur 
if Im[pij] > 0. 

In the next section we present approximate solutions of Eq.s ( |13a| )-( |T3l|) , both the temporal and 
the steady state cases. These will be compared with numerical calculations of the full system, 
i.e., without any approximation. 

IV. DENSITY MATRIX EQUATIONS IN THE DRESSED STATE BASIS IN THE 

SECULAR APPROXIMATION 

As mentioned before the Hamiltonian part of the master equation has a simple form in the 
dressed state basis given by Eqs. ( paH9c|) (the Hamiltonian is diagonalized in the dressed state 
representation). The problem arises when the spontaneous emission and pump terms are present 
in the master equation,Eq.(|l|), giving the complicated couplings appearing in equations ( |13a| ) - 



3i|). Solving exactly the complete set seems to be a formidable task even with "Mathematica" 
|26| . However, the situation can be simplified if the frequency difference between the dressed 
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states of the manifold, namely the Rabi flopping frequency R is large compared with the rates 
7 b , 7 C ,A. We can then ignore the "nonsecular" terms, i.e., couplings between populations and 



coherences (see [f25[|). 

a. The evolution of the population terms. When the "nonsecular" couplings between popu- 
lations and coherences are ignored we obtain from ( |13a| ), ( |13d| ) and ( |131| ) the following equations 
for the populations: 

Paa = -(r Q + A')p aa + -A'ppp + ^A'p 77 , (15a) 

ppp = -(Tfs + -A)pp P + -{T a + A') Paa + (I> + \a- -A> 77 , (15b) 

Pll = -( F /3 + ^ A )P77 + ^(r Q + A')Paa + (I> + ^A - ~A')ppp. (15c) 



Note that population conservation is still maintained. The interpretation of equations ( |15a|) - 
( |15c| ) is very clear. The state \a > is depopulated with a rate (T a + A'), which in turn, is 
distributed equally among the states \(3 > and (7 >, as can be readily seen from the one 
half factor multiplying the coefficient of p aa in ( |15b| ) and ( |15c| ). The state \a > is also being 
populated with a rate ^A by the states \j3 > and (7 >. The state \(3 > (I7 >) is depopulated 
at a rate (Tp + |A) and repopulated with the same rate from (7 > (|/3 >). The set ( |15a| ) -( |15c| ) 



can be solved exactly by calculating its eigenvalues and eigenstates, subject to the condition 
Paa + P/3/3 + P 77 = 1- This yields the temporal solution 

Paa(t) = Pta + \Pa*(0) ~ P^]^^ (16a) 
PPfi® = P% + \PW{0) - P% + \(Paa(0) ~ pl)]e- (2rd+A ^ A,)i 

- \\Paa{0) -P^e-^^, (16b) 

Prrit) = P 77 ~ ~ P% + \(PU0) ~ p s L)]e~ {2Tri+K -^' ]t 

-\[p a a^)-p s L]e- {V ^ h,)t . (16c) 

where p«(0) are the initial populations. The steady state populations are given by 

A' 

p ™ = 2r Q + 3A" (17a) 
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We can see that the population ppp, p 77 have similar behavior, though not identical, a 
fact which is not surprising at all in light of the very similar composition of the states \/3 > 
and 1 7 >. The population of the state \a > is unique in the sense that it decays with only 
one decay constant, while the other populations have a composite decay Figure 3 shows 
a comparison between the exact solution for p aa and ppp (solid line), obtained by solving 
numerically the equation set (|13a|) -( |13f|) , with our approximate analytic solution (|16a|) - (|16cj ) 



shown in dashed line. The normalized parameters for the numerical simulation were set to be 
Q = 207c, 7b = 27c, G = 0-ljc, A = 37 c . We can see that the population p aa is a monotonically 
increasing, oscillating function of time which reaches a steady state value p s ^ a ~ 0.27. The 
behavior of ppp is opposite, i.e., it is a monotonically decreasing oscillating function and it 
reaches the steady state value p% ~ 0.36. p 77 is not shown because its similarity to ppp, 
due the choice of parameters made. The approximate solutions describe nicely the envelope of 
oscillation and the correct expression for the steady state. One can see that p aa < ppp = p 77 , 
for any finite A, and hence population inversion do exist in the dressed state basis. For the 
transitions \a >— > \/3 >, \a >— > (7 >, \/3 >— > \j > and (7 >— > |/3 > the population difference 
is negative (noninversion) and it remains to be seen whether this transitions amplify and thus 
result in lasing without inversion in the dressed state basis. In the strong coupling field limit, 
Q >> G the steady state population become 

7c + A 

^-^r-2^T3A' (18) 

and 

Paa = 2^T3A- (19) 

In the following we will get into more detail regarding gain without inversion in the dressed 
state basis. 

b. Evolution of coherences. Ignoring the "non secular" couplings between coherences and 
populations in Eqs. (|13a|) - (|13fl ) results in equations that are simpler than the original ones, 
however, they are still very complicated, particularly the equations for /3g 7 and it's conjugate, 
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which are coupled to all the other coherences. Hence one would like to further approximate 
these equations in such a way that the resulting solutions will be fairly simple on one hand, and 
be a reasonable approximation to the exact solution on the other. We solved numerically the 
complete set (|13a| ) - (|13t| ) and found that p^ and hence its conjugate are substantially larger 



than the other coherences, indicating the crucial role these coherences play. In light of the 

above, we couple each atomic coherence to itself (describing the free evolution) and to p/3 7 and 
p 7/ 3 acting as the dominant source terms. 
This gives the following equations: 

Pa/3 = -(r Q/3 - iR)p a /3 ~ f'Pfa + (f - F)p 7 /3, (20a) 

p ay = -(r Q7 + iR)p ay + (f - f')p/3 7 - f> 7/ 3, (20b) 

p^ = -(I> 7 + 2iR) Pfh - (Tfs + 1-A- \m)p iP , (20c) 



along with the equation for p 7/3 = p*^ . Solving the eigenvalue problem of eq.'s ( |20a| ) - (|20c| ) 
we find the transient solutions for the coherences , in the strong coupling field limit. These 
solutions are: 

x n „ (T -f')(T + A/2- A'/2) + i4f> R 

= Aexp[ - (r ^ - ' R)t] + c ( r, + A/2-A72)(r, 7 -r e , g + M) e ^~^ + M « 
+ J (f--r )+; P ( r a + A/2 A 7 2)/4 fl ^ _ 2im (2la) 

f '{Yb + A/2 - A72) - iAR(V - F) 
p Q7 = £exph(r> + + - A/2 _ A//2)(r ^ _ ^ - ^ expRF^ + 2iR)t] 

f' -i(T -r)(T g + A/2- A72)/4R r „. . , 

+ 1 r - p ' ' 11 exp[-(r^ 7 - 2xJ2)f], (21b) 

p = c< 4 ^ - exp[-(r« 7 + 2iR)t\ + D i Tp + A/2 ~ AV2 x 

^ 7 I> + A/2-A'/2 PL v ^ Ji AR 

X exp[-(r Bj -2iR)t], (21c) 

where A, B, C and D are constants ought to be calculated from initial conditions. Figure 
4 shows the exact coherences obtained by solving numerically eq.'s (|13a|) - (|13f ) (solid line), 



and the corresponding approximate solutions based on the analytic expressions (|21a|) - (|21cj) . 
The chosen parameters are the same as in Fig 3. It can be seen that the exact solutions reach 
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a steady state different from zero, while our approximate analytic solution have zero steady 
state values. This is due to the absence of the population source terms in the approximation. 
The oscillation frequency is predicted correctly by the approximate solutions. Moreover, the 
approximate solution for p^ appears to be more accurate than the other two, again indicating 
the crucial role played by p^ and p 7j g. Note also the negligible contribution of p a p and p ai to 
P/3 7 . The latter coherence is not coupled to p a p, p ai yet the approximation remains satisfactory. 
Eq. (|21cj) also shows that p^ has an almost pure sinusoidal form of frequency 2R. To the 
contrary p a /3,Pai3 and p ai have a composite oscillation, being a superposition of frequencies. It 
can be seen from the numerical solution presented in figure 4 that the coherences p a /3, and p ai 
possess a definite sign, thus the corresponding transitions between dressed states, being either 
amplified or attenuated. More precisely Im{p a p) is positive in all the range shown, hence the 
transition \(3 >^ \at > is amplified (with population inversion in the dressed state picture at the 
frequencies VLl + R, fl p + R ). The transition \a >— > (7 > is also amplified since Im(p ia ) > 0, 
however, without population inversion in this case. This situation is very different from that 
occurring in the bare state basis, where the coherences oscillate back and forth across zero, thus 
experiencing periodic amplification and absorption see Fig ^, and [T2j . In contrast to p a p and 
/Vp the sign of Im{p 1 p) is alternating thus the transition \(3 >—> \j > is being amplified and 
absorbed periodically. Another feature reminiscent of figure 4 is the strength of the coherence 
P/3 7 , which is seen to be three orders of magnitude stronger than the other two coherences. The 
most striking result is that of the transition (7 >— > \a >. It is absorbing despite the population 
inversion (see figure 4 (b) ). 

The main deficiency of eq.'s (|21a ) - ( picj ) is the zero steady state predicted by them. The 
reason for this is the omission of the source terms (populations) in writing eq.'s ( |20a|) -( |20cj) . 
We have solved analytically Eqs. (|20a]) - (|20cj ) with the source terms included. The solution 
obtained was checked against numerical calculations and found to be in excellent agreement. 
Unfortunately, the solution is so complicated, that even reduction operations carried out by 
"Mathematica" |2B[ could not give a manageable solution. For the purpose of finding the steady 
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state coherences it is sufficient to retain the population terms in ( |20a| ) -( |20c|) , set to zero the 
time derivatives, and solve the resulting algebraic equations. This yield the following steady 
state coherences: 

St _ I St \* _ + ^)[T a P + (£q + 2A')(f + f')] 

Paf3 {Paj) (r2 /3 + J R2)( 2 r Q + 3A') 

ATp{T a p + iR){Y a + AQ[(2F - Y){2Y P + A - A'/2) - 2MT] 
{Y% + R*)(2Y a + 3^)^ + 4^-^ + A/2 -A'/2)«] ' 1 &J 

st 4r^(r a + AQ (2T> + A - A 7 / 2) - 2»fl 

P/37 2r a + 3A' (r /3 + A/2-A'/2) 2 -(r2 7 + 4 J R 2 )- 1 ' 

The steady state values predicted by the last expressions were checked against numerical cal- 
culation and found to be in very good agreement. Note that expressions (|22a|) - (|22b|) also 



give the steady state dispersion and not only the gain or absorption coefficients. The expected 
dominance of p^ on the other two coherences, can be seen by noting that the term r + Y', 
appearing in p a p and p Q7 (but not in p^) varies like o(G/R), thus the coherence p a p varies 
like p a /3 ~ G/R 2 while p^ varies like p^ ~ 1/R. The second term in Im{p a p) is of the order 
G/R? and can be neglected. Amplification of the \(3 >— > \a > transition at frequencies Ql + R, 
fl p + R, occurs whenever Im(p a p) > 0. Taking into account only the first term in ( |22a| ) we find 
gain in the following two cases: 

1. For any incoherent pump rate A (even zero), if > 7c- The physical interpretation of 
this result is clear: if level \b > is drained more quickly than |c > there is no need for the 
incoherent pump, and the recycling of the population is accomplished by the coherent 
probe field. 

2. For n ^ R2 7e < lb < 7c provided that the incoherent pump rate is strong enough such that 



A > 



r Q (7e-7b) 



r Q -2Q2( 7c _ 7f) )/R2- 

This gain is "regular" gain, due to population inversion since ppp > p aa . The transition 
\a >— > 1 7 > at frequencies Ql + R, Q p + R will be amplified for the same range of parameters 
because the corresponding gain coefficient is proportional to Im(p 7a ) = —Im{p ai ) = Im(p a p). 
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However this transition is inversionless, and it is due to external field induced quantum inter- 
ferences and atomic coherences. The opposite transition, I7 >— > \a > is at frequencies Ql — R, 
flp — R. It exhibits absorption with population inversion (see figure 4 (b)). In a sense this is the 
reverse process of amplification without inversion and it is explained as a constructive quantum 
interference for the stimulated absorption process. The \(3 >^ I7 > transition, at frequencies 
Ql + 2-R, fl p + 2R will be absorbed, since the term 4R 2 appearing in the denominator of (|22b| ) 
far exceeds the other denominator terms , resulting in Im{p^p) < 0, and hence absorption. 
The transition I7 >^ \f3 > at frequencies Ql — 2-R, fl p — 2R in turn will be amplified for any 
incoherent pump rate. 

In the strong field limit the imaginary parts of the steady state coherences can be expressed 
in terms of the original atomic parameters as 



Wapj ^yaJ 2 v / 2fi 2 (2 7c + 3A) V ; 

MO - i 7 '!!l (23b) 



As mentioned before, in the strong coupling field limit (week probe) the state \a > becomes 
the highest excited bare state |c >, thus we are facing a situation of full noninversion in the 
dressed state basis, i.e., p cc < ppp = p 77 (see Eq.. ([18|) - (|T9|)). Eq ( |23a| ) indicates that gain can 
be obtained for the |c >— > \(3 > and |c >— > I7 > transitions (the Autler- Townes transitions 



Tl\ ) for the following conditions: 

1. For any incoherent pump rate A if 7& > 7 C . 

2. For ±7 C < 7b < 7 C , provided that A > > (7c ~ 7 '' ) . 

This gain is without inversion. The physics involved in this condition is as follows: 
It states that the dephasing time of level \c > namely (|7 C ) _1 must be longer than any other 
decay time in the system. The dephasing process must be slow with respect to other processes 
in order to preserve the phase of the dipole transition p ac . This makes possible the quantum 
coherent effect whereby the interference can result in gain without inversion. 
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To improve the temporal results for the populations we retained only the dominant coher- 
ences, namely the approximate solutions for p^ 7 and p 7/ g in the population equations (Eq. (|15a|) - 
(|15c|) serving as source terms. These equations are integrated giving the following: 



-(T +§A')t 



p 77 (t) = p 



st 
77 



+ c 2 e 



(2r /3 +A-iA')i _ ±nP ar 



-c\e 



-(r 0+ |A')t _ C2e -(2i>+A-iA0t _ 



1 

2> 



(24a) 
(24b) 
(24c) 



where the particular solution is given by: 

pZ(t) - 



P/3 7 (0)A' (T a + §A' - T^) cos(2Rt) + 2Rsm(2Rt) 



(r Q + |A'-r /?7 ) 2 + 4i? 2 



e _r /37* 



The integration constants are given in terms of initial populations and coherences by: 

cx = Paa (0) - - pZ(0) 

c 2 = Pw(0) - p% + ^(p aa (0) - p s * a ) 

Figure || shows the difference between the exact population p aa and the approximate solution 
( p4a|) - (|24cj) . It can be seen that the approximate solution is very accurate. 

To compare with the steady state situation in the bare state basis, we need to transform 
back to the bare state basis ,via the the matrix product p B = T~ 1 p Dr T, where p Dr is the density 
matrix in the dressed state basis, formed by the steady state populations and coherences of 
Eq..'s ( |17a| ) -( |17b|) , and ( |22a|) -( |22b|) . Further utilization of the strong coupling field limit gives 
bare state populations, in the steady state regime. 

lc + A 



Paa 

and 

Pec 



Pbb 



7c + 3A : 



A 

7c + 3A' 



(25a) 



(25b) 



The bare state coherences are: 
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lb (7c + A) 



(26a) 



Pab 



—I 



2fi(2r c + 3A)' 
Y A( 7& - 7e) - 7c 
r 2fi 2 (2 7c + 3A) ' 



Pac 



%G 



(26b) 



and 



G 7c 



(26c) 



Pbc 



Q(2 7c + 3A) ' 



In obtaining the expressions for populations and coherences we find that our general form 
reduce to previously obtained results |i2| . 

We see that in the bare state basis one always has 7m(p a &) < 0, thus the coupling laser is 
always attenuated. The probe transition exhibits inversionless gain for 75 > 7c for pump rates 



true lasing without population inversion can be realized, both in the bare state and the dressed 
state basis. 



Absorption in the presence of inversion and amplification without inversion in a three level 
V - type system are found in the dressed state picture. Both of these effects are the manifes- 
tation of the quantum interference that occurs in multilevel systems. Moreover, the above two 
processes constitute a manifestation of a complementarity principle. 

We have presented an analysis of light amplification without population inversion in this 
system within the framework of the dressed state basis. The equations of motion for the 
elements of the density matrix are derived from the master equation. In the dressed state 
picture new relaxation terms are defined that are related dierctly to the coherently prepared 
states and the quantum interfernce effects, interferference terms are identified. They are shown 
to be the source for amplification without inversion and for absorption despite the inversion. 
Consequently, approximate analytical time dependent solutions for dressed state populations 
and coherences were obtained. Comparison of these approximate solutions with the numerically 
calculated quantities shows excellent agreement. Both of these solutions exhibit the familiar 



satisfying A > 



7i>-7> 




. From the analysis presented above we conclude that for a week probe , 



V. SUMMARY 
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Rabi oscillations. Steady state density matrix elements were also calculated, from which we have 
concluded that for a weak probe field, true lasing without inversion exist for the appropriate 
incoherent pump rate, i.e., lasing without inversion in any state basis. Steady state quantities 
were transformed back to the bare state basis, and were found to be in perfect agreement with 
results in the literature. Conditions for inversionless gain where obtained and show physically 
sound basis. Finally, the new feature of absorption despite population inversion found in this 
calculation emphasizes the importance of quantum interference. Quantum interference is thus 
shown to imprint its effect on the processes in the dressed state picture. 
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|(3(N + 1)(N'+1)> 

|a(N+i)(N'+i)> 

|Y(N + 1)(N'+1)) 
|(3(N)(N'+1)> 

|a(N)(N'+i)> 
Iy(N)(N'+d) 

I (3(N - 1)(N* +1) > 
la(N-l)(N'+l)> 
I y(N -1)(N'+1)> 



|P(N + 1)(N')> 

|(X(N + 1)(N')) 
|y(N + l)(N')> 

IP(N)(N')> 

|a(N)(N')> 

lY(N)(N')) 



I (3(N - 1)(N' ) > 
I a(N - 1)(N' ) > 
I y(N - 1)(N' ) > 



FIG. 2. Manifolds e(N,N'),e(N+l,N'),e(N,N' + 1), and e(N + 1, N' + 1) etc. of uncoupled states 
of the atom + lasers photons (left hand part). The dressed level (perturbed levels) are shown at the 
right hand side. The circles represent steady state populations. 
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FIG. 3. Normalized time evolution numerical simulation, of dressed state population p aa and pgg 
obtained by solving Eq.s (13a) - ( |13f| ) (solid line), and the approximate solution based on Eq.s ( |16a| ) 
-(|16cj) (dashed line) .The chosen parameters are: f2 = 207 c , 7j, = 2^ c , G = 0.l7 c , A = 37 c 
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FIG. 4. Normalized time evolution numerical simulation, of dressed state coherences obtained by 
solving Eq.'s ( |13a| ) - (|13f ) (solid line), and the approximate solution based on Eqs. ( |21a| ) - (|21c[ ) (dashed 
line). The chosen parameters are the same as in Fig || 

. Note the absorption despite population inversion seen in (b). 
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FIG. 5. Difference between Exact population p aa obtained by solving Eqs. ( |l3a[) - ( |l3f|) and the 
approximate solutions based on solutions (24a-[24c), 
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